Some Considerations on the Additional Absorption Peak in the c-axis Infrared 
Conductivity of Bilayer Cuprate Superconductors: Interpretation of the Changes 
Induced by a Parallel Magnetic Field and the Role of Bilayer Splitting 
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Changes of the 400 cm~^ peak in the c-axis conductivity of underdoped YBa2Cu306.6 upon ap- 
phcation of a parallel magnetic field reported by Kojima et al. are shown to be consistent with the 
model where the peak is due to the superfluid. Results of our calculations of the c-axis response 
of bilayer compounds with well defined bilayer-split bands are presented and discussed. For mod- 
erate values of the bilayer splitting (Ae comparable to 2Amax) the spectra of the superconducting 
state exhibit an additional mode which is due to the condensate and similar to the one of earlier 
phenonomenological approaches. 



The c-axis infrared spectra of bilayer high-Tc cuprate 
superconductors (HTCS), i.e., HTC compounds with 
two Cu02 planes within a unit cell, exhibit the follow- 
ing interesting features. At low temperatures a broad 
absorption peak develops in the spectra of the real 
part of the optical conductivity cti in the frequency re- 
gion between 300 cm~^ and 800 cm~^ (its frequency is 
material- and doping-dependent)i*2i^*^. At the same 
time some of the phonon modes are strongly renor- 
malized, i.e., their frequencies and/or linewidths and/or 
spectral weights changeiSiSi^. The magnitude of the ad- 
ditional peak and the changes of the phonons are par- 
ticularly spectacular for the strongly underdoped super- 
conductor YBa2Cu307_5 (Y-123)3i&. For a long time 
these phenomena had retained their mystery. In 2000 
Griininger et al^ showed that the additional peak in 
the spectra of Y-123 can be understood and quantita- 
tively described in terms of a model proposed in 1996 
by Van der Marel and TsvetkoAjS. Within this model, 
the so called multilayer model, a bilayer compound is 
considered as a superlattice of homogeneously charged 
superconducting planes and interbilayer and intrabilayer 
regions as shown in Fig. 1. 

The model dielectric function £{lo) is given by 
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where £oo is the interband dielectric constant and the 
conductivities au and Oint are defined by ju = cfuEu 
and jint — cTintEint- In the superconducting state (SCS), 
they contain contributions of the condensate. 



hi j int 



(3) 



where lou and LOint are the plasma frequencies. Further 
and u^^f are the regular parts of the conductivi- 
ties. In the absence of the latter, eiuj) exhibits two zero 
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FIG. I: Schematic representation of the multilayer model. 



crossings, at wli = \/uiint /soo and at ujl2 = \/ww/eoo, 
corresponding to two longitudinal plasma modes: the in- 
terbilayer and the intrabilayer one; in addition, it has a 
pole at 



{dbi + dint)Soo 



(4) 



corresponding to a so called transverse plasmon (TP). Its 
spectral weight is given by 



St — -^£0 



{du + dint){dblU)fnt + dint^li) 
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Griininger et al. attributed the additional peak to the 
TP. The model has then been extended^ by including 
local fields acting on the ions. This allowed the authors 
of Ref. to explain not only the additional peak but 
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also the related phonon anomalies. Today, the multilayer 
model^ and the local field idea™ are becoming accepted. 
The nature of the intra-bilayer currents responsible for 
the additional peak, however, is the subject of an ongoing 
discussion^'^. Here we address two related issues. 

Changes of the additional peak in strongly underdoped 
Y-123 induced by a parallel magnetic fieldh&Ye been stud- 
ied by Kojima et al}^ . It has been found that the field 
suppresses the condensate at w = 0. Its spectral weight 
(SW) is redistributed between a so called a mode at ca 
40cm~^ and the additional peak at ca 400 cm~^. The 
SW increase of the latter amounts to ca 460f7~^cm~^. 
Surprisingly, the field does not influence the peak fre- 
quency. Below we show that this behaviour is consistent 
with the claim that the peak corresponds to the super- 
conducting condensate. Assume for simplicity that the 
suppression of the superfluid in the interbilayer regions 
is complete, i.e., ojint — 0, that the field does not mod- 
ify the intrabilayer plasma frequency uju ( uju » ^int), 
and that the regular parts of the two conductivities can 
be neglected. The field induced changes of St and uut can 
then be obtained using Eq. (5) and Eq. (4), respectively. 
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It appears that St somewhat increases and ojt very 
slightly decreases. For the values of the parameters 
obtained by fitting the spectra of a similarly under- 
doped Y-123 sample, uu = 1200 cm^^, uoint = 220 cm^^ 
(Ref. 0), du = 3.3 A and d^t = 8.4 A, we obtain 
ASt = 800fi"^cm"2 ^nd Aojt = -3cm~^ In reality, 
the suppression of the condensate in the interbilayer re- 
gions is certainly incomplete. Kojima and coworkers suc- 
cessfully interpreted the low-frequency part of the spectra 
in terms of a model where the Josephson coupling is con- 
siderably weakened only in ca 50% of the interbilayer re- 
gions. Using this model we obtain ASt ~ 400ri^^cm^^ 
and Alut ~ — 1.5cm~^ [the magnitudes are reduced ca 
by a factor of 0.5 with respect to the simple estimates 
of Eqs. (6) and (7)}ii. The value of ASt is in excellent 
agreement with experiment. The very small value of Aujt 
accounts for the fact that, at first glance, the 400 cm^^ 
peak does not shift upon application of the field. 

Transverse plasmon in the presence of bilayer split- 
ting. The multilayer model of Eq. (1) is well justifiedi^ 
for a system of weakly ("Josephson") coupled Cu02 
planes where t±^ max is much less than any important in- 
plane energy scale, in particular t_Lmax < Amax- Here 

max and Aniax are the maximum values of the inter- 
planar hopping matrix element and the superconduct- 
ing gap, respectively. The overdoped bilayer compound 
Bi2Sr2CaCu208 (Bi-2212), however, exhibits two bands: 
the bonding band (BB) and the antibonding band (AB), 
with the maximum splitting of ca 90 meV correspond- 
ing to t_Lmax — 45 mcV > Amajj^^*^- Here we present 



very preliminary results of our calculations of the c-axis 
response in such a system. Details will be presented in 
Ref. Il5i The main ingredients of our approach are sum- 
marized below. 

Our starting point is the tig ht-binding formalism out- 
lined in Appendix B of Ref. [la (with a more general form 
of the interplanar hopping matrix elements). For sim- 
plicity, the spacing layers will be supposed to be fully 
insulating (i.e., t±i„i = in the notation of Ref. H^ . 
First, we express the intrabilayer conductivity uuiy) de- 
fined by jui'jj) = (^bi{<jj)EM{Lo)- It is given by Eq. (B9) of 
Ref. 
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where a is the in-plane lattice constant, T is the intrabi- 
layer c-axis kinetic energy per unit cell. 



{T) = - 



2a^ 



dk t±_ (k) (riBfe - UAk ) , (9) 
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and X is the relevant current-current correlation function. 



x(q,w) 



di([j(q,t),j(-q,o)])e(t)e^' 



(10) 

The interplanar hopping matrix element t±^ is related to 
the dispersion relations of the two bands as 2tj^(k) = 
EAfe ~ ^Bk ■ The occupation factors are denoted by UBk 
and UAk- In Eq. (10), j(q, t) is the Fourier component of 
the paramagnetic part of the intrabilayer current density 
operator. It can be shown^^ that 



j(q = 0) 



^ t^{]f){c+j^CAK - c\jiCBK) , 



K=(k,fc^) 



(11) 

where c^^, cbk^ c\j^, and cak are the conventional 
quasiparticle operators. The Matsubara counterpart 
xiiuJm) of x(q = 0,w) is given hyi^ 



46^4; 1 
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J dk4(k)x 



X < -^^Tr[GB(k,ip„)GA(k,ip„ + iL^„0] } , (12) 

where Gb and Ga are the Nambu propagators. Note 
that for Gb = Ga (i.e., in the limit of weak interlayer 
coupling) we recover Eq. (10) of Ref. m When deriving 
Eq. (12), vertex corrections have been neglected. In or- 
der to express the dielectric function, we have to consider 
interlayer Coulomb interactions. Random phase approx- 
imation leads to Eqs. (1), (2) with aint = 0. The mul- 
tilayer formula (1) is thus justified even if the coupling 
between the closely-spaced planes is strong (the coupling 
across the spacing layer, however, must be negligible). 
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The results presented below have been obtained us- 
ing the above equations with selfenergy corrections ne- 
glected. The final formula for x reads 



e^dh 



dk<i(k)) 
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huj + i6 — E, 



hu> + id + E, 



where 



^i/2(k) 



^Bk^Ak + ^Bk^Ak T EskEAk 



(13) 



(14) 



EBkEAk 

^Bk/Ak are the superconducting gaps, Egk/Ak = 
\J^'Bk/Ak + ^sfc/Afc' Fermi function and Ef = 

EAk i EBk- A conventional form of the dispersion 
relationsiii has been used: 

eB/A{k) — ~2t[cos{kxa) + cos{kya)] — 



At' cos{kxa) cos{kya) =F 



t±, 



[cos(fc2;a) — cos(fcj,a)]^ — ii , 
(15) 

and the bands have been assumed to possess the same 
superconducting gap 



A 



Bk — ^Ak 



[cos(fca;a) — cos(fcj,a)] . (16) 



Figure 2 shows the main results 

for the following values of the model parameters: 
dhi = 3.4 A, dint — 12.0 A (the values correspond to 
Bi-2212), Eoo = 5.0, t = 250 meV, t' = -100 mcV, 
*_Lmax = 45me\4i^, fi — — 350meV, Amax = 30mcV, 
and S = 1.0 meV. In the normal state (NS, Amax — 
0, T = 100 K), £bi2{^) contains only the interband 
(IB) contribution. Its maximum is located at lujb ~ 
2t_Lmax/?i. In the SCS (T = 20 K), ebi2iuj) consists of 
the IB (pair-breaking) contribution centered at lujb ~ 
^V^imax + ^max/^> with a Smaller SW than in the NS, 
and the one of the condensate (C) at zero frequency. 
The ratio SW(C)/SW(IB) depends on Amax/i_Lmax- For 
Amax = 20meV (30meV; 60meV), e.g., we obtain 
Ubi = 1870 cm-\ SW(C)/SW(IB) = 0.47 (2260 cm"!, 
0.93; 2690 cm^i, 3.07). The conductivity cri(w) has 
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FIG. 2: Calculated spectra of the real (ebii) and imaginary 
(£612) parts of the "intrabilayer dielectric function" eti defined 
by Eq. (2) (a) and of the real part of the c-axis conductivity 
(b). 



sharp maxima labelled as P and Q (we do not con- 
sider here the small maximum labelled as Q'). They 
occur at frequencies up and ujq, for which etii = 
— {du / dint)£ 00 and ebi2 is small [see Eq. (1)]. The max- 
imum Q has a large SW of about 60000 f)~^cm~^. It 
appears already in the NS and it will be broadened 
by selfenergy effects. The large difference between loq 
and LOiB is caused by the interlayer Coulomb interac- 
tions. The maximum (mode) P appears only in the 
SCS as a result of broken gauge symmetry. Its fre- 
quency and SW {Sp) increase dramatically with increas- 
ing Amax/iimax. For Amax = 20meV (30meV; 60meV), 
e.g., we obtain ujp — 350 cm~^, Sp 
(490cm~i, 250017-icm-2; 840 cm^^. 
For Amax/ij_max ^ (i.e., in the limit of strong inter- 
layer coupling), the mode vanishes. Obviously, the maxi- 
mum P can be related to the TP of the phenomenological 
model involving a superlattice of interbilayer and intrabi- 
layer Josephson junctions. The underlying excitation can 
be visualized as a resonant oscillation of the condensate 
density between the closely-spaced Cu02 planes. 
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In conclusion, the changes of the 400 cm"'^ peak in- 
duced by a parallel mg. field reported by Kojima et al. 
are consistent with the claim that it is related to the su- 
perconducting condensate. The regular parts of the two 
conductivities need not be taken into account in order to 
explain the changes. The c-axis response of bilayer com- 
pounds has been studied considering the B-A splitting. 
The conductivity cti exhibits a maximum in MIR related 
to the B-A transitions. Its shape will be determined by 



selfenergy terms and it can be expected to change slightly 
when going from the NS to the SCS. Only in the SCS, an 
additional maximum appears in FIR, related to the con- 
densate, with similar properties as the one of the naive 
Josephson-superlattice model. Our results indicate that 
the latter model is compatible with moderate values of 
the B-A splitting. 
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1 C. C. Homes et al, Physica C 254, 265 (1995). 

2 J. Sclmtzman et al., Phys. Rev. B 52, 13665 (1995). 
^ C. Bernhard et al, Phys. Rev. B 61, 618 (2000). 

* V. Zelezny et al., Phys. Rev. B 63, 060502 (2001). 

^ S. V. Dordevic et al, Phys. Rev. B 69, 094511 (2004). 
® D. Munzar et al. Solid State Commun. 112, 365 (1999). 
M. Griininger et al, Phys. Rev. Lett. 84, 1575 (2000). 

* D. van der Marel and A. Tsvetkov, Czech. J. Phys. 46, 
3165 (1996); Phys. Rev. B 64, 024530 (2001). 

® T. Timusk and C. C. Homes, Solid State Commun. 126, 
63 (2003). 



^° K. M. Kojima et al, Phys. Rev. Lett. 89, 247001 (2002). 

J. Marek and D. Munzar, to be published. 
12 N. Shah and A. J. Millis, Physical Review B 65, 024506 

(2001). 

1^ D. L. Feng it et al., Phys. Rev. Lett. 86, 5550 (2001). 
" Y. D. Chuang et al, Phys. Rev. Lett. 87, 117002 (2001). 

1® D. Munzar et al, to be published. 
1® D. Munzar et al, Phys. Rev. B 64, 024523 (2001). 
" O. K. Andersen et al, J. Phys. Chem. Solids 56, 1573 
(1995). 



